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ABSTRACT 

We present a first-principles theory of Resonant Relaxation (RR) of a low 
mass stellar system o rbiting a more massive black hole (MBH). We first extend 
the kinetic theory of iGilberfj (|l968l) to include the Keplerian field of a black 
hole of mass M,. Specializing to a Kepl erian stellar system of mass M <C M ,, 
we use the orbit-averaging method of ISridhar fc Toiimal (j2ni5: Paper ll ) to 
derive a kinetic equation for RR. This describes the collisional evolution of a 
system of A ^ 1 Gaussian Rings in a reduced 5-dim space, under the com¬ 
bined actions of self-gravity, 1 PN and 1.5 PN relativistic effects of the MBH 
and an arbitrary external potential. In general geometries RR is driven by 
both apsidal and nodal resonances, so the distinction between scalar-RR and 
vector-RR disappears. The system passes through a sequence of quasi-steady 
secular collisionless equilibria, driven by irreversible 2-Ring correlations that 
accrue through gravitational interactions, both direct and collective. This cor¬ 
relation function is related to a ‘wake function’, which is the linear response of 
the system to the perturbation of a chosen Ring. The wake function is easier 
to appreciate, and satisfies a simpler equation, than the correlation function. 
We discuss general implications for the interplay of secular dynamics and non¬ 
equilibrium statistical mechanics in the evolution of Keplerian stellar systems 
toward secular thermodynamic equilibria, and set the stage for applications 
to the RR of axisymmetric discs in Paper III. 


Key words: galaxies: kinematics and dynamics—galaxies: nuclei 


1 INTRODUCTION 


Star clusters around massive black holes (MBH) in galactic nuclei are dense stellar systems in which 
stellar orbital motions are fast enough for relaxation processes to have significantly modified their 
structure. Relaxation is driven by gravitational “collisions” between stars, leading to exchanges of 
orbital energies and angular momenta. Classical two-body relaxation results in the diffusion of both 
energy and an gular momentum, over th e two-body relaxation time scale Treiax (IChandrasekhar 


I942I . Il943al jbl: iBinnev &: Tremaind l2008l ). In its fully developed form two-body relaxation can 


feed the MBH with stars, and lead to the formation of a stellar density cusp around the MBH 
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2 Sridhar & Touma 


(jBahcall &: WolfI Il976l : ICohn &: KulsrudI Il978l ). But for the stellar systems observed in galactic 
nuclei, Trpiay can ex c eed th e Hubble time; then two-body relaxation would have had negligible effect. 
Rauch Sz Tremain^ ( 1996l l — hereafter RT96 — proposed a more efficient mechanism, Resonant 
Relaxation (RR), that could be in operation in stellar systems with globally degenerate (or resonant) 
orbital frequencies. RR does not enhance the energy relaxation rate but promotes more efficient 
relaxation of orbital angular momentum distribution. When this happens in a low mass stellar 
system orbiting a MBH, stars whose orbits have become sufficiently eccentric can fall into the 
MBH, feeding it mass, energy and angular momentum. RR is particularly important in a low-mass 
stellar system around a MBH, and is a result of the degeneracy of the Kepler problem. As the radial 
and azimuthal periods of a Keplerian ellipse are equal, orbital precession is suppressed; only the 
smaller non-Keplerian forces, such as cluster self-gravity, general relativistic effects and external 
gravitational sources contribute to orbital precession. This gives rise to a pattern of interactions 
between Keplerian orbits which persists over the longer time scales of precession. It is this extended 
period of interplay, or coherence, between Keplerian orbits that makes RR more efficient than 
classical two-body relaxation. 

Sketch of the Rauch—Tremaine model: Let us consider a stellar system of size R and total mass 
M, consisting oi N ^ 1 stars (each of mass m* = M/N), orbiting a MBH of mass M, . Stellar orbits 
are determined by the combined gravitational effects of the Keplerian potential of the MBH and 
cluster self-gravity, general relativistic effects and any external gravitat ional sources. The stellar 
system will be referred to as Keplerian (jSridhar &: Toumal 12015: Paper iT l when the non-Keplerian 
forces are much smaller than the Keplerian held of the MBH. Therefore, for an isolated Keplerian 
stellar system the mass ratio must be small, e = M/M, <C 1, and its size must be large enough, 
R ^ r, , where r, = 2GM,/(? is the Schwarzschild radius of the MBH. In two-body relaxation 
random gravitational encounters between stars, with a coherence time of order the short Kepler 
orbital timescale T^ep = 27r (ii^/GM,)^^^ , build up to order unity over times Treiax = 0 

As noted earlier T^eiax can exceed the Hubble time, whereas RR works more efficiently because 
of the longer coherence times which we now estimate. Following RT96 we consider RR in a non- 
relativistic Keplerian stellar system, for which [r,/R) <C e <C 1. Then orbits evolve slowly over 
secular timescales T^ec = £~^7kep j which is longer than Kepler orbital times. Over several Tsec, 
stellar orbits can be thought of as Gaussian Rings (see Paper I), which are Keplerian ellipses of 
fixed semi-major axes that can deform in the mean self-gravitational field of the cluster. A secular 
collisionless theory of this system was worked out in Paper I, and applies in the limit N ^ oo while 
M = Aim* is held constant. For realistic systems N is large but finite, say N 10® — 10^, and it 
is this granularity that drives RR. 

As the semi-major axes of Gaussian Rings do not change, there is no exchange of Keplerian ener¬ 
gies during secular gravitational encounters □ But mutual torquing can exchange angular momentum 
between Rings leading to diffusion of angular momenta. The diffusion of a vectorial quantity like 
angular momentum can be complicated, and depends on the geometry and orbital structure of the 
stellar system. According to RT06 the physics of this process can be understood by considering two 
extreme cases: scalar-RR which is the diffusion of the magnitude of angular momenta (and hence 
the eccentricity of stellar orbits); and vector-RR which is the diffusion of the direction of angular 
momenta. Diffusion of the magnitude of angular momenta occurs as a result of the mutual torquing 
of Gaussian Rings whose apsides precess over times Tgec- Hence the coherence time for these secular 
encounters is ~ Tseo which is much longer than T^ep- The magnitude of the angular momentum 


^ Logarithmic corrections to T^eiax have been neglected. 

^ The Keplerian orbital energies diffuse over the longer time scale T^eiax. 
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can be thought of as varying linearly in time for times shorter than Tgec, and random-walking over 
longer times. The timescale for scalar-RR turns out to be = HT^ec = ^T^elax j which is shorter 
than Ti-eiax- Vector-RR can occur in near-spherical systems, even in the absence of a MBH, when 
the apsides of stellar orbits precess fast (say, when the stellar system is immersed in an external den¬ 
sity cusp) and apsidal resonances have negligible effect. Then orbital eccentricities are conserved, in 
addition to the semi-major axes. Averaging over the fast apsidal precession, each stellar orbit can be 
replaced by an axisymmetric annulus whose inner/outer radii are equal to the (constant) peri/apo 
centre distances. Diffusion of the orientation of orbits occurs as a result of the mutual torquing of 
the stellar annuli with coherence times, ~ In a Keplerian stellar system it turns out that 

the vector-RR timescale is also of order the coherence time, giving = N^/‘^Tsec = 
which is shorter than even the scalar-RR timescale. Thus there emerges the notion of a hierarchy 
of relaxation processes, with longer coherence times corresponding to shorter relaxation times: (a) 
the Keplerian orbital timescale and classical two-body relaxation; (b) the slow secular apse preces¬ 
sion timescale and scalar-RR; (c) the slower secular nodal precession timescale and vector-RR. For 
general stellar distributions the distinction between scalar-RR and vector-RR can get blurred. 


Resonant Relaxation over the years: RR has been explored extensively through both nu¬ 
merical simulations and stochastic modeling, a trend that was also initiated by RT96: restricted 
A^-Body simulations (splitting stars into field and test stars) and A^-wire simulations were used to 
demonstrate that RR is more efficient than two-body relaxation; the efficiency was characterized in 
both coherent/lin ear and diffusive regimes, by htting for parameters that were left free in their ran¬ 
dom walk model. iHopman &: Alexanden (j2006l l revisited the problem of the feeding of a MBH with 
RR acting as the driv ing mechanism, instead of the classical two-body relaxation considered by 
Bahcall fc WolfI ( 19761 ). Using the RT96 model they parametrized the strength of RR in a Fokker- 


Planck equation (which is averaged over angular momentum) for the distribution of energies, with 
loss cone boundary conditions. They could constrain the fluxes into the MBH, the energies of stars 
that contribute to it, and enhancement in gravitational wave events. They also studied the implica¬ 
tions of the process for th e Galactic Centre and conc luded that the kinematics of stars was consistent 
with the workings of RR. iGiirkan fc HopmanI (120071 ) used the wire-approximation of RT96 and ran¬ 
dom sampling to explore the RR ti me scale as a function of eccentr icity, and concluded that RR 
is more efficient on eccentric orbits. iKupi. Eilon fc Alexander! (|2010l l used small-scale Newtonian 
A^-body simulations to determine the efficiency of RR with regard to the inspiral of compact rem¬ 
nants into the MBH and emission of gravitational waves, and concluded that RR will increase the 
rat es of inspiral events by a fac t or of a few over t wo-body relaxation. 


Madigan. Hopman fc Levin ( 2011 1 generalized Hopman fc Alexander ( 2006l l in two significant 


ways by using: (a) An autoregressive moving average model for RR (instead of using the RT06 
random walk model) that is constrained by restricted A^-body simulations; (b) A Fokker-Planck 
equation that kept track of both energy and angular momentum. Their solution allowed them to 
account for a density core in the Galactic Center, and to conclude that binar y disruption makes 
for re sonantly relaxed orbits which are too eccentric for the observed S-stars. iKocsis fc Tremaine 
(j2nill ) explored an analytical model of vector-RR to investigate possib le warps in the distr i bution 
of stars at the Galactic Centre. Vector-RR has also been studied by iKocsis fc Tremaind (j2015l ) 
through numerical simulations of annuli-annuli interactions, and a random-walk model of orbit 
normals on a sphere. They found that the dominant torques were between stars with radially over- 
lapping orbits, and t hat th e vector-RR rate increased rapidly for highly eccentric (e > 0.8) orbits. 
Touma fc Tremaine ( 20141 ) determined maximum entropy equilibria of self-gravitating Keplerian 


discs, for the case where all the particles have equal semi-major axes. Their numerical investigations 
revealed broken-symmetry states which could be lopsided and uniformly processing. 
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4 Sridhar & Touma 


Using extensive post-Newtonian N^-body simulations. iMerritt et al.l (|201lh . argued that general 
relativistic precession offers a severe barrier to the capture of stars by the MBH, thereby quenching 
RR. They considered dynamical mechanisms for barrier penetration and assessed the rate of gravi¬ 
tational wave events they enable. This is the most serious attempt at studying the workings of RR 
through A^-body simulations, and stands as a benchmark for co nfrontation with any funda, r nenta. l 
theory of RR in the presence of general relativistic corrections. iBar-Or &: Alexander] (j2014l . l2015h 
have also studied the problem of the quenching of RR torques by general relativistic precession; they 
model the background potential as a correlated Gaussian nois e, and study the stochastic process 
throug h a Fokker-Planck equation and Monte Carlo simulations. Hamers. Portegies Zwart fc Merritt! 
(j20l4 l made extensive post-Newtonian, restricted N^-body simulations with a view to understand¬ 
ing the behavior of S-stars in the Galactic Centre. They identified three regimes of relaxation: 
non-resonant, resonant and anomalous, and recovered functional forms for the diffusion coefficients 
in all three regimes. These results were used in a Fokker-Pl anck equation to obtain the steady-state 
distribution of angular momentum for orbits near the MBH. Merritt ( 2ni5a b 3l developed a numer¬ 
ical algorithm of the Cohn-Kulsrud type for solving a Fokker-Planck equation in both energy and 
angular momentum, which used diffusion coefficients describing the effects of two-body relaxation, 
RR, anomalous relaxation and energy loss due to emission of gravitational waves; loss of stars to the 
MBH was also included. The algorithm was applied to recover time-evolving distribution functions 
in galactic nuclei. 

Approach to a Kinetic Theory of Resonant Relaxation: The brief overview of developments 
in RR given above is, admittedly, far from exhaustive; rather it is meant to be illustrative of what 
we perceive as some major trends over the past two decades. Numerical simulations have grown in 
sophistication; in quality and size, as well as in their ability to include physical effects such as general 
relativistic precession and emission of gravitational waves. In contrast theoretical understanding of 
RR is based on various parametrized stochastic models whose parameters are recovered from the 
numerical simulations, ft is necessary to develop a new framework for RR, one whose foundations 
are based on rigorous theories of non-equilibrium statistical mechanics. We begin by recalling the 
theoretical basis underlying the standard treatment of classical two-body relaxation. 

Two-body relaxation has been studied extensively both analytically and numerically. The 
historical beginnings of the k inetic theory of stella r system s are rooted in the problem of the 
evolution of globular clusters. Chandr asekhad (j 19421 . Il943al jbl) derived a Fokker-Planck equation 
by considering gravitational encounters in an infinite homogene ous system of stars. The deriva ¬ 
tion of the diffusion and frictio n coeffi c ients was simplified by Cohe n. Soitze r fc Routlv ( 1950l h 
Gasiorowicz. Neuman fc Riddell (1956); Rosenbluth. MacDonald fc Judd ( 1957rL The general ki¬ 
netic equation for inhomogeneous stellar systems was derived by iGilbertI (1 19681 ): this is rigor¬ 
ous and based on a systematic 1/N expansion of the BBGKY ( Bogoliubov , Born , Green, Kirk¬ 
wood, Yvon) equations of non-equilibrium statistical mechanics. Hevvaerts ( 201Cll ) revisited the 
problem using action-angle variables (assuming an integrable system), and proved that the ki¬ 
netic equation satisfies an H-theorem. The most difficult parts of the theory are certain polariza¬ 
tion terms arising from collective inte r action s. A first simplification of the kinetic equation was 
achi eved by iPolvachenko &: ShukhmanI (jl982l i through the neglect of the polarization terms — see 
also Chavanis ( 2012 . 2013 ll^ There has been some work on the applications of the kinetic equa¬ 
tion (in the Gilbert form, or the reduced Polyachenko-Shukhman form), but we do not review 
this literature. We note, however, that the kinetic equation applies to inhomogeneous systems for 


® This simplification is analogous to the passage from the equations of 
equations, for electrostatic plasmas. 


Balescu 




and 


Lenard 


( igeoll 


to the 


Landau 


( igfidi 
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Stellar Dynamics around a Massive Black Hole II 5 


which computing interactions bet ween realistic stellar orbits can be complicated. One can pass from 
Polvachenko &: ShukhmanI (|1982l ) to the more tractable and well-investigated Fokker-Planck de¬ 
scription of classical two-body relaxation. This involves the step taken originally by Chandrasekhar, 
in which gravitational encounters are estimated by assuming that stars move with constant veloci¬ 
ties. Thu s we are comforted that th e journeyman theory of two-bod y relaxation, w hich is textbook 


via 


ma terial (IBinnev &: Tremaine 


Polvachenko h ShukhmanI (|l982h . 


20081) ■ rests on the firm foundations of iGilbertI (|l968D through descent 


But we have seen that two-body relaxation is not efficient in galactic nuclei with MBHs, and we 
must consider RR. To the best of our knowledge, the stochastic models used to interpret numerical 
si mulations of RR are all ad hoc. Our goal is to develop a first-principles theory of RR, that is based 


on 


GilbertI (j 19681 ). For simplicity we consider a stellar system consisting of equal mass stars; the 


extension to a stellar system with a range of masses is straightforward, but this will not be treated 
here. Our focus is on presenting the derivation of the kinetic equation governing RR in detail, so all 
the developments in this paper are of a formal nature. Applications of our theory will be presented 
in separate papers. We begin in § 2 with an outline of Gilbert’s theory of the collisional relaxation 
of general stellar systems, adding to it arbitrary external sources of gravity. The theory is based on 
a 1/N expansion: at order unity we get the collisionless Boltzmann equation (CBE), and collision 
terms emerge at order 1 /N. We cast the equations in Poisson Bracket form by defining various 
potentials. In § 3 the external potential is taken to be that of a MBH, and a canonical transformation 
is effected to coordinates centered on the MBH. Another canonical transformation expresses the 
Gilbert equations in the Delaunay variables, which are the natural action-angle variables for the 
Kepler problem. No approximations have been made and the equations apply to a stellar system 
around a MBH for any mass ratio e. In § 4 we consider the application to Keplerian stellar systems 
for which e <C 1. We begin with the collisionless theory and recall the results of Paper I on the 
secular limit, obtained using the method of multiple scales. Then we move to the collisional theory 
and orbit-average the Gilbert equations; the calculations are straightforward, with details given in 
the Appendix. The result is a kinetic equation for RR as envisioned by RT96. In § 5 we supplement 
the mean-self gravity of the system by secular relativistic corrections (up to 1.5 post-Newtonian 
order) and external sources of gravity. Thus we arrive at a general kinetic equation for RR. Some 
general physical features of this equation are discussed in § 6, and concluding remarks are offered 
7. 


m 


2 GENERAL SETTING FOR COLLISIONAL RELAXATION 
2.1 Outline of Gilbert’s theory 


We begin with a brief account of the theory of collisional relaxation of [Gilbert] (Il968l ). The stellar 
system (or cluster) consists of a large number N of stars, each of mass m*, and total mass M = Nm^,. 
Let Xi and Vi be the position and velocity of the fth star (for i = 1,..., A), with respect to an 
inertial frame. The phase space of the system is 6A-dim with coordinates Fj = {xi,Vi). We also 
use dFj = d?Xid?Vi and write the 6A-dim volume element as dri...dr 7 v- We recall that, in 
the standard statistical formulation (iLifshitz fc Pitaevskiilll98ll ). the system is described by an N- 
particle distribution function (DF) (Fi,..., Fjv, t) which is a symmetric function of Fi,..., F^r. 
Reduced s-particle DFs ... ,Ts,t) derived for s = 1,..., (A—1), by integrating over 

the coordinates Fg+i,..., Fjv. The A-particle DF is a probability distribution in the 6A-dim phase 
space {Fi,..., Fat}; hence all the s-particle DFs are probability distributions in their respective 
6s-dim phase spaces {Fi,... ,r 5 }. Evolution equations for the /(*) can be obtained, by integrating 
the Liouville equation over these very coordinates. This is the BBGKY hierarchy of equations, in 
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6 Sridhar & Touma 


which the equation for depends on f^‘^\ and so on up to For a large N system the chain 

of equations is very long, and is usually truncated through ordering in a physically relevant small 

parameter. __ __ 

Drawing on earlier work in plasma physics (|Balesculll960l : lLenardlll960f ). Gilbert achieved clo¬ 
sure of the BBGKY equations for stellar systems through power-series expansion in the small 
parameter (1/A^). With notation F = {x,v), F' = {x',v') and F" = {x'',v"), the 2-particle DF 
can be expanded as: /(^)(F,r',t) = , where 

is the irreducible part of the 2-particle correlation, satisfying f dF = f /W dF' = 0. It is 
straightforward to verify that is determined uniquely in terms of and to first order 
in 1/N. Hence the BBGKY hierarchy closes at this order, and provides two coupled equations for 
and which describe the collisional evolution of a self-gravitating system. Irreversibility 
arises through the “adiabatic turn-on” initial conditions, for which /W vanishes in the distant past. 
Henceforth we drop the subscript “1”, and write the 1-particle DF as f{T,t). The acceleration of a 
test particle at {x,t) by the mean-field is: 


a{x,t) = / f{r',t)b{x,x')dT', where h{x,x') = GM- 


X — X 


\x' — x\ 


( 1 ) 


is the (scaled) Newtonian “bare” inter-particle acceleration. External sources of gravity (a MBH, 
an in-fall ing gas c loud or star cluster) can also be included in the formalism, although these are 
absent in [Gilbert] (1 19681 ). Let a^{x,t) be the acceleration of the test star due to external sources. 
Then the equation for f(x,v,t) can be written as: 


dt 


+ + (a + a®) 

ox 


dv 


1 df 


= —a 


N dv 




dv 


dr'. 


( 2 ) 


Mean-field theory (MFT) corresponds to the limiting case when Y” —>■ oo, m* —?■ 0 with Ym* = 
M = mass of the system which is held constant. The accelerations b(x,x'), a(x,t) and a^(x,t) 
are all well-defined in this limit. However, both terms on the right hand side vanish in the limit. 
Hence the MFT of a stellar system is completely described by the 1-particle DF, f{x,v,t), whose 
time evolution is governed by the collisionless Boltzmann equation (CBE). When Y is large but 
not infinite, the accrual of small changes induced in /, by the right hand side of eqn.([2]), cannot 
be ignored over long times. The first term corrects the mean-field acceleration of a star by a small 
amount, from a to (1 — 1/Y)a, taking note of the fact that the gravitational force on a star is due 

only to the other (Y — 1) stars. It is the second term that is of importance for collisional evolution, 
( 2 ) 

where acts as the driving term. 

The equation for {T,T',t) depends on both f{T,t) and f(T',t), and is difficult to deal with 

in the general case. However, we are interested in describing the quasi-steady collisional relaxation 
of /. The stellar system can be thought of as evolving through a series of collisionless equilibria, 
for which {v-d/dx + ((1 — 1/Y)a + a®) -d/dv}/ ~ 0, being of order 1/N or smaller. This implies 
that df/dt is hrst order in 1/Y, due solely to the collision term. It turns out that, in this case, 
can be expressed in terms of a conditional probability function ^(r | r',t), which will henceforth 
be referred to as the Wake of F' at F. The wake can be dehned through the following gedanken 
experiment: from the Y stars that are distributed in 6-dim phase space according to /, select one 
star and place it at the phase space location F'; this will induce a small perturbation (1/Y)VV at 
every location F in the entire phase space. At this point it is useful to define two kinds of acceleration 
at any point x. The wake acceleration. 


a'^{x,T',t) 


j b{x,x'')W{T'' \ T\t) dF", 


(3) 
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is the acceleration at x due to the wake of F'. We also define the perturbing acceleration, 

aP{x,x\t) = b{x,x') — a{x,t), 


(4) 


where the mean-field acceleration is subtracted from the bare inter-particle acceleration. The 
wake ^(rlr'jt) must be calculated by integrating an equation for the related wake function 
W(r I r'(T), T), where r'(T) = {x'{t')) is the location at time t' ^ t, of the star that ar¬ 
rives at the desired location T' at time t. The equation for W(r | is: 

+ V-— - h {a{x,t) + a {x,t)} )~ 


dt' 


dx 


dv dv 

= —a'^{x,x'{t'),t')'^, fort' ^ t. 


(5) 


Eqn.([5]) needs to be solved with the adiabatic turn-on initial condition VV(r | r'(T), t') —)• 0 as 
t' — oo, corresponding to a vanishingly small wake in the distant past. The right hand side 
identifies the perturbing acceleration a^{x, x'{t'),t') as the driver of the wake. The wake VV(r | T', t) 
can be obtained by (formally) integrating eqn.Q from t' = —oo to t' = t, and is therefore a 
functional of the quasi-steady /. It may be verified that eqn.Q preserves J" W(r | T', t) dT = 0, 
because the wake contains no net mass. 

Extending the wo r k of 


Rostokeii (j 19641 1 in plasma physics to the case of an inhomogeneous 


stellar system. iGilbertI (|l968h used non-trivial operator identities to derive the fundamental relation 
( 2 ) 

between /. and W: 


= w(r|r',t)/(r',t) + w(r'|r,t)/(r,t) 

+ [ >v(r I r", t) w(r' | r", t) /(r", t) dr". 


( 6 ) 


(‘ 2 '\ 

This decomposition of means that the irreducible 2-particle correlation at (T, T') gets three 
contributions from the wake function: (a) The wake of T' at T; (b) The wake of T at T'; (c) 
The product of the wake values at the points T and T' of a third Ring at T", summed over all 
locations IZ". All three contributions come with suitable /-weighting. The third term accounts 
for the contribution of collective effects (“gravitational polarization”) to the microscopic processes 
driving RR. This is the nature of the full theory at 0(1/A^). Since f T', t)dr' = 0, eqn.Q 

implies that f ^(r | T', t) /(R, t)dr' = 0. 

Together, eqns.Q and ([6]), give a closed set of kinetic equations governing the collisional relax¬ 
ation of /(r,t). This can also be written more explicitly as: 

I + + {(' - v) “+-I = 


c^'V] = -^Ib{x,x').-^{f{r,tmr'\r,t)} dr'. 


(7b) 


C^'^V] = -^ I /(^^ t) {b{x, x') + a-(a:, R, t)} • dR . (7c) 

Here and are the dissipative and fluctuating parts of the collision term. As the 

notation indicates both are functionals of /, because they depend on the wake W(r | T',t) which, 
as discussed earlier, can be thought of as a functional of /. 
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8 Sridhar & Touma 


2.2 Gilbert’s equations in Poisson Bracket form 

Gilbert’s kinetic theory is expressed in the physically direct language of the various forces acting 
on a star. However, it is more convenient to exchange accelerations for corresponding potentials, 
and rewrite the kinetic equations in Poisson Bracket (PB) form. We have, 

{b, a, a*, aP} = - A{p, n, H", HP} , (8) 

where the potential functions on the right hand side are defined by: 


P{x, x') 

GM 

1 /I ’ 

\x — x'\ 

(scaled) Inter-particle Poisson kernel; 

(9a) 

n(ic, t) 

= I fiT',t)P{x,x' 

)dr'. 

Mean field potential of the cluster; 

(9b) 

U%x,t) 

= External potential; 


(9c) 

u-{x,r',t) 

II 

|r',t) dP' 

' , Wake potential at x due to P'; 

(9d) 

nP(a;, x', t) 

II 

1 

t), 

Perturbing potential. 

(9e) 

The mean-field Hamiltonian is: 





= Y + n(a;,t) + U%x,t), (10) 

with the corresponding mean-field equations of motion. 


dtc 

dt 




= V, 


dt> 

dt 




= a + a . 


dv dt dx 

The PB between two phase space functions xi(r,t) and X 2 (r,t) is: 

r , def dxi <9x2 dxi dx2 

Xl ) X 2 (r) g^ 

Using this we can rewrite eqn.([5]) for the wake function ^(r | T'{t'),t') as: 

+ [/(r,t'),n-(a^,r'(t'),0](r) 

= [np(®,®'(t'),t')>/(r,t')](r) > ^ 


dt' 


+ 


with adiabatic turn-on initial condition lim ^(r | = 0. 

t'^ — OO 


( 11 ) 


( 12 ) 


(13) 


There are two ways to write the Gilbert kinetic equations in PB form. The first is based on 
eqns.dJ]) and ([6]), which is closer to the original BBGKY equations: 


dt 


+ 


/(r,t),F«(r,t)- 


n(a;, t) 
N 


(r) 






/i 7 (r,r',t) = >v(r|r',t)/(r',t) + >v(r'|r,t)/(r,t) 

+ y'>v(r|r",t)>v(r'|r",t)/(r",t)dr". (u) 

The second is the explicit form given by Gilbert, displaying the dissipative and fluctuating contri- 
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Stellar Dynamics around a Massive Black Hole II 9 


butions to the collision integral. From eqns. (f7a|) - (l7^ we have: 

= , 


dt 




(r) 


1 


I [P{x,x'),f{T,t)W{r'\T,t)]^^^ dC, 


(r) 


(15a) 

(15b) 


c^"V] = ^//(r',t)[P(*,®') + n"'(a;,r',i),>v(r|r',i)](^)dr'. (i5c) 

Either form of the kinetic equations applies to the collisional evolution of / whenever it satisfies 

~ 0, being of order 1/N or smaller. (16) 


/(r,t),i/W(r,i)-II(M 


J(r) 


Therefore /(r,t) may be thought of as passing through a sequence of collisionless equilibria in a 
quasi-static manner. 


3 KINETIC EQUATION FOR A STELLAR SYSTEM ORBITING A MASSIVE 
BLACK HOLE 

We now specialize the kinetic theory to a stellar system of total mass M = Nnii,, orbiting a MBH 
of mass M,. The mass ratio e = M/M, is a natural ordering parameter; the smaller it is the more 
“Keplerian” is the stellar system. We begin with the kinetic equation in the BBGKY form, eqn. (ll4jl . 
and make two canonical transformations. The first transformation is from the inertial frame to a 
(generally non-inertial) frame centered on the MBH. Then we transform to the Delaunay variables, 
which are natural action-angle variables for the Kepler problem. The formulae we derive are valid 
for arbitrary e; however the second transformation proves really useful only when e is small. 


3.1 Transformation to coordinates centered on the MBH 


The Keplerian force of the MBH on a star is described by taking the external potential of ean. (l9^ 
as: 


U%x,t) = 


GM, 

X - x,{t)\ ’ 


(17) 


where x,{t) is the position of the MBH with respect to the inertial frame. It is natural to transform 
to a non-inertial frame centered on the MBH. Let r = x — x,{t) and u = v — x,{t) be the 
position and velocity, respectively, of a test star relative to the MBH. It is evident that {r,u) 
are a canonically conjugate pair of 6-dim phase space coordinates. Using notation T* = {r,u) 
we have dT = d^aid^t) = d^rd^rr = dT*. We write the (1-particle) DF as /(r*,t) and the 
wake as W(r* | r(., t). We take as an 0(1) reference the Kepler potential of the MBH potential 
jqext _ given in ean. ffTTll . Then the other potentials in eqns.([9a|), (l9b)) . (I9dp and (f9e)) are 
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10 Sridhar & Touma 


all quantities of 0(e). Hence we rescale these and write {P, n,n",nP} = e{p, (/j™, where 
GM, 


p{r,r') 

p{r,t) 


\r — r' 


, (rescaled) Poisson kernel; 

/(r(^, t)p{r, r') dr(,, (rescaled) Mean field potential; 


= jp{ry)WK\T'„tWi, (rescaled) Wake potenfal; 

p^{r,r',t) = p{r,r') — (p(r,t), (rescaled) Perturbing potential. 


(18a) 

(18b) 

(18c) 

(18d) 


All the Poisson Brackets in ean. (ll3p and eans. (ll5all - (ll5cll remain invariant, because {x,v) —>■ 
(r, u) is a canonical transformation, i.e. for any two phase space functions xi and X 2 y the PB defined 
earlier in ean. (ll2p is: 


r 1 _ r 1 def dxi dx2 dxi dx2 

[xi,X2](r) - [xi,X2](ro - 


(19) 


du dr 

The first term on the left hand side of the wake eqn. fll.Sp is {dW/dt')(^r) where we have now indicated 
explicitly that P = (x, v) is to be held constant. We need to express this in terms of functions of 
P* = {r,u) and t. The chain rule gives, 


dt' ) (p) 


W ) 


— X 


(r*) 


dW 

* dr 


— X 


.. ^ 
* du 


V dt' ) (P^) 


W , r-x,{t) - u-x,{t) - -|i,(f)p 


J{r*) 


( 20 ) 


Then ean. ffTsp for the wake function W(r* | is: 

^ + [w,Porg(r*,0](r.) + e[/(r*,0,<^"'(or'(t'),0](r.) 

= e [c/?P(r,r'(t'),t') 1 /(r*,t')](r,) > ^ 


with adiabatic turn-on initial condition lim VV(r* | r*(t^), t^) = 0. 

t'^—oo 


( 21 ) 


We have written Porg(^’) ■*.1) t) = Hd)-\-{r-x,{t) — u-x,{t) — \x,{t)\‘^/2}, and dropped the subscript 
P* in the time derivative. The acceleration of the MBH, x,{t ), is due to the gravitational forces of 
all the stars, and is a 0{e) quantity. Rescaling, we write x,(t) = eA,{t), where 


A,it) = GM. I /(r*,t)^dr*. (22) 

Using eqn. (flop for Hd) we can write 

u? GM 

Porg(r*,f) = Horgir,u,t) = — - - + e(p{r,t) + er-A,{t) (23) 

2 r 

explicitly as a function of {r,u,t). This is the Hamiltonian governing the mean-field dynamics of 
r and u through the equations of motion: 


dr 

dt 


dH, 


org 


du 


= u. 


drr 

dt 


dH, 


org 


GM, 


dr 


-r — £ 


dp 

dr 


eA. , 


(24) 


as can be verified from the definitions of r and u. Note that both Porg and the equations of motion 
are identical to those given in § 2 of Paper I. Similarly, ean. dldp can be rewritten in the MBH 
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Stellar Dynamics around a Massive Black Hole II 11 


variables. Then /(r*,t) satisfies: 


dt 


/(r*,t), i^org(r*,t) - —v?(r,t) 

■/V J (r*) 


N 




J(r*) 


= >v(r*|r:„t)/(r',t) + >v(n|r*,t)/(r*,t) 

+ [ >v(r* I T':, t) >v(n I r:, t) /(r:, t) dr 


(25) 


3.2 Transformation to the Delaunay variables 


We now transform from the MBH-centric canonical coordinates, T* = {r,u), to the Delaunay 
variables P = {J, L, fej , which are action-angle va riables for the Kepler problem (jPlummer 

1960 : Murray Sz Dermott 1999 : Binnev Sz Tremaine 2008 1. The three actions are: I = \/GM,a ; 


L = I\/l — the magnitude of the angular momentum; and = Lcosi the z-component 
of the angular momentum. The angles conjugate to them are, respectively: w the orbital phase; 
g the argument of the periapse; and h the longitude of the ascending node. The transformation 
is exact and no restriction is placed on the magnitude of e (although this is preparatory to the 
orbit-averaging that follows in the next section for e <C 1). In the Hamiltonian Hoj-g of ean. (|23p . 
the Keplerian orbital energy E]^ = {u^12 — GM,/r^ will be taken as 0(1). The other two terms, 
Eip and er-A,, are evidently of 0(e). 

The Kepler orbital energy Ei^{I) = —l/2(OM,//)^ depends only on the action I. Therefore, 
for the pure Kepler problem (i.e. when e = 0) all the Delaunay variables, excepting w, are constant 
in time; w itself advances at the rate 


Dk(/) 


dl 


{GM,f 

P 


(26) 


In order to rewrite the various physical quantities in terms of the T) variables we need to express 


the position vector rj'D'] = jx, y , z) in terms of the T> variables (IPlummeiill96fll : iMurrav fc Dermott 


19991 : ISambhus fc Sridhaii 120001 ): 

/ X \ / GgGh - G^ShSg -SgGh " GiS^Gg S,Sh \ / a{Gg - e) \ 

y = GgSh + GiGhSg -SgSh + G,GhGg -SiGh aVT^Sg 

\ z j Y SiSg SiGg Gi j \ 0 y 


(27) 


where S and G are shorthand for sine and cosine of the angles given as subscript. Here a is the 
semi-major axis; g is the eccentric anomaly, related to the orbital phase through w = {g — eshig ); 
e = Vl - is the eccentricity; and i is the inclination angle determined by cosf = {Lz/L). It 

is also useful to note that r = = a(l — ecosr/). 

The DF is written as f(T>,t) and the wake as W(P | D',t). The transformation T* —)• P being 
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12 Sridhar & Touma 


canonical, dF* = dP. Then the potentials of ean. (|18a11 - (ll8dh and A, can be written as: 
GM, 


p{V,V') 

A,{t) 

pP{V,V',t) 


\r — r' 


(rescaled) Poisson kernel; 


f{'D',t)p{T>,T>')dT)', (rescaled) Mean field potential; 


GM, J f{'D,t) — d'D, (rescaled) MBH acceleration; 

J p(P, 'D")yV{'D" I V', t) dV" , (rescaled) Wake potential; 

p{T>,T>') — , (rescaled) Perturbing potential. 


The mean-field Hamiltonian is: 

Horg{V,t) = Pk(-f) + eip{V,t) + er-A,{t) 


(28a) 

(28b) 

(28c) 

(28d) 

(28e) 

(29) 


Since the transformation T* —)• P is canonical, all the PBs remain invariant. The PB between any 
two phase space functions Xi(^) X 2 (P) is: 

def f dxi dx2 dxi <9x2 A f dxi dx2 dxi <9x2 


[Xl, X2](6) ( 


dw dl dl dw 


dg dL dL dg 


dh dLy dLy dh 


(30) 


From eon. (1211) . we see that the wake function yV{'D \ satisfies 


BW BW 

+ 9" 1 d- r-A, + e [ /(P, t'), ip^{V,V'{t'),t') ] 


Bt' 


( 6 ) 


= s[(pP{V,V'{t'),t') , , forp ^ t, 


with adiabatic turn-on initial condition lim yV{'D\'D'{t'),t') = 0. 
Similarly, from eon. (1251) . the kinetic equation for /(P,t) is: 


(31) 


^ n ^ 

m + " 


/> ( 1-]^ ) + r-A, 


N 


1 ( 6 ) 


p{v,v'), 


1 ( 6 ) 


dP'. 


/l//(P,P',t) = W(P|P',f)/(P',t) + W(P'|P,f)/(P,f) 

+ [ >V(P I P", t) W{V' I P", t) /(P", t) dP" , 


(32) 


4 ORBIT-AVERAGING FOR A KEPLERIAN STELLAR SYSTEM 

The stellar system is said to be Keplerian when the mass ratio e = M/M, <C 1. In this case the 
dominant force on a star is the inverse-squared Newtonian force of the MBH, and e is the natural 
small parameter for a perturbation theory of Keplerian systems. We want to simplify the description 
of the previous section and obtain an 0(e) description of the quasi-steady collisional evolution. 
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Stellar Dynamics around a Massive Black Hole II 13 


4.1 Elements of Secular Collisionless dynamics 


The first step is to understand the general features of the collisionless limit, N ^ oo and m* ^ 0 
with M = constant, of the kinetic ean. (l32p for f(T>,t): 

^ + e[/> </^](6) = 0- (33) 

We begin by recalling results from the theory presented in Paper I, which is accurate to 0(e). 
The natural time scales in the problem are (a) the Kepler orbital period Tkep = on which 

the phase w varies, and (b) the longer secular time scale Tgec = e~^Tkep over which the other 5 
Delaunay variables 7^ = {/, L, L,-,g, h} evolve in the mean0 The CBE describing secular evolution 
was derived by using the method of multiple scales, to orbit-average the full CBE over the fast 
orbital phase w. Two time variables t and a slow time t = et were introduced; these are natural 
measures of time over T^ep and Tgec; respectively. We argued that the secular description cannot 
depend on the fast orbital phase w. The remaining 5 Delaunay variables IZ = {I,LjI'z]g,h} are 
the coordinates of a Gaussian Ring (or just Ring) which is a Keplerian ellipse of given semi-major 
axis, eccentricity, inclination, periapse angle and nodal longitude. The secular evolution ol N ^ 1 
stars can be thought of as the evolution of ^ 1 Gaussian Rings, so the 0(1) part of the DF must 
be a function of the slow variables and independent of the fast variables {w,t). It turns out 

that the appropriate expansion of the DF, to 0(e) accuracy, is: 

^ {F(7^,r)+e/o(7^,r)} + ^ ^ ^(7^, t, r) exp [inu;] + 0(6^), (34) 

n^O 


where F{IZ, r) will be referred to as the Ring DF. The DF / is a probability distribution function in 
P-space satisfying f f dV = 1 in the absence of loss of stars to the MBH — see eqn.(l) of Paper I. 
Using dT> = dIZdw , we have 


F(7^,r)d7^ + e J /o(7^,r)d7^ + 0{e^) = 1. 


(35) 


The method of multiple scales provided a solution to 0(e). We returned to the usual single-time 
description, and presented the following solution to the full CBE eqn. (j33p : 

f{V,t) = ^{F(7^,e^)+e/o(7^,e^)} + ^AF{V,et) + w - n^t, et) + 0(8^). (36) 

ZTT ZTT ZTT 

Here F{IZ, r) satishes the Ring CBE: 


dF 


(37) 


where <h(7^,r) is the Ring Hamiltonian (coming from stellar self-gravity), 

dm 

^ ’ 

which is equal to the orbit-averaged mean-held potential. The PB without subscript is the 4-dim 
PB, 


^{IZ,t) = y (/?(T’,t) 


(38) 


[Xi, X 2 ] = 


dxi 9x2 dxi 9x2 


+ 


d)ad)^ _ 9xi 9x2 

dh dLy dL:, dh 


(39) 


dg dL dL dg 

whose action is restricted to the 4-dim I = constant surfaces in the 5-dim 7^-space. Under the 
self-consistent Ring dynamics of ean. (l37p . F(IZ,t) can have signihcant time variations over times 


There are three kinds of 0(e) corrections to the Ring description. fo{IZ, r) is a secular correction 


^ The relevant time scales for the two best—known Keplerian systems, the Galactic Centre and the nucleus of M31, are discussed 
at the end of the introductory part of § 2 in Paper L 
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term to the Ring DF, which remains undetermined in our 0(e) theory. The other two correction 
terms are purely fluctuating in w. Of these, AF{'D,£t) is a known function of the Ring DF F and 
w; every Ring DF F{TZ, et) in the 5-dim 7^-space is accompanied by a small fluctuating (in w) 
distortion in the full 6-dim P-space. A is an arbitrary function of its arguments TZ, {w — 
and et, with zero mean over w. The arbitrariness of form is due to the limits of our 0(e) theory. 
Yet it has physical significance, in that it emphasizes the requirement that there are no instabilities 
growing on the fast orbital times T^ep — as it must be, for F(TZ, r) to be a physical 0(1) description. 

The 0(1) self-consistent secular collisionless description is obtained by taking the limit e ^ 0, 
while the time scales of interest are large enough to make r a sensible 0(1) measure of time. This 
limit is described completely by the Ring DF, F(TZ,t), which satisfies the Ring CBE (I37p . The 
normalization of ean. (f35ll reduces to 

J F(n,T)dn = 1, (40) 


so F(TZ,t) may be thought of as a probability distribution function in the 5-dim 77-space. The 
orbit of every Ring is restricted to its 4-dim I = constant surface. The dynamics of the remaining 
4 Ring variables {L, L^', g, h} is given by the Hamiltonian equations, 

dL (9<k d^ 54* dL^ 94> dh 94> 
dg ’ 


(41) 


dr dg ’ dr 9L ’ dr dh ’ dr dL^ ' 

which will be referred to as the Ring equations of motion. Using these, it is straightforward to write 
the the Ring CBE ((371) as dF/dr = 0, justifying the term “collisionless”. The Ring description 
makes no reference to the orbital phase w, and is self-consistent and closed in 5-dim 77-space. 
However, we will later need to take account of the evolution of w to 0(1). It turns out that this is 
slaved to the Ring degrees of freedom, and given by: 


drc 

dr 


Dk(/) 


+ 


( 94 > 


(42) 


The orbital phase increases steadily over times Tkep; a-nd also has modulations over times Tg; 


4.2 Orbit—averaged kinetic equation 

We now return to the problem of quasi-steady collisional evolution of Keplerian stellar systems, 
and simplify the kinetic equation of the previous section through first order perturbation theory in 
e. In quasi-steady evolution it is clear that the DE / cannot depend on t, and hence must be of 
the form f('D,T). Moreover, in the collisionless limit, the DE must go smoothly to the form given 
in eon. (1341) . Therefore we write: 

- OO 

^ X] /n(^,T)exp [irew] O(e^), (43) 

n=—oo 

where F ^ 0 and /_„ = f* . We recall from § 2 that the wake is a “response function” characterizing 
the collisionless system which, to 0(1), is described completely by the secular DE F{JZ,t). Hence 
to 0(1) the wake must describe the perturbation in 77-space due to the selection of a Ring at a 
given location. Therefore the power-series in e for the wake can be written as: 

. OO 

yV(V\V' ,t) = —kF(77|77',r) -|- ^ ro„(77 | P', r) exp [inrc] -|- 0(e^), (44) 

ZTT ZTT 

n=—OO 

where {\/N)W (77 | 77', r) can be thought of as the additional perturbation at 77, resulting from the 
selection of a Ring at 77'. We will refer to W (77 | 77', r) as the Ring Wake function. 

The Gilbert equations, eans. (|3ip and (|32p . are valid for any value of the mass ratio e. Orbit¬ 
averaging is a systematic, perturbative method that can be applied to a Keplerian stellar system 
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for which e <C 1. This proceeds by substituting the multiple-scale forms, ean. (l43p for the DF / 
and ean. (|44l) for the wake W, in the Gilbert equations and developing them perturbatively in the 
small paramater e. The calculations are straightforward yet tedious, and we provide the details in 
the Appendix. Henceforth we discuss the 0(1) theory for the collisional evolution of F(IZ,t) and 
W(IZ I IZ', t) in 5-dim 7^-space — the fluctuating quantities, fn and Wn, contribute only at higher 
order. In contrast to Gilbert’s kinetic theory of point masses this is a kinetic theory of Gaussian 
Rings, and is the appropriate framework for describing the Resonant Relaxation of RT96. Analogous 
to the various potentials that arise in the PB form of Gilbert’s equations, various Ring potentials 
make their appearance. These are all derived in the Appendix — see eans. llAlTbl) . (lAlSh . ()A20h 
and (IA2ip — and reproduced below: 


nn.n') = 


1 


r — r 


is the “bare” inter-Ring potential in 7^-space, and plays a role analogous to the Poisson kernel, 
p{'D,'D'), in Gilbert’s theory. The other Ring potential are: 


4>(7^,T) 


J F{IZ',t) 'I'(7^, 7^')d7^', Ring mean-held potential; 




J W{IZ" I 7^^, r) T {IZ, IZ") dIZ" , Ring wake potential; 


4>P(7^,7^^T) 


'I/{IZ,IZ') — ^{IZ,t) , Ring perturbing potential. 


All of these have the same signihcance as the corresponding potentials in Gilbert’s theory, as is 
clear from their appearance in the governing equations for W and F given below. 

From eqn. (IA22p the wake function kF(7^ | 7^'(r'), r^) is determined by: 

dW r nr 

— + [W{IZ\IZ'{t'),t'),^IZ,t')] + [F(7^,T'),^"(7^,7^'(r'),r')] 


= [4>P(7^,7^'(T'),r'), T(7^,r')] , for r' ^ r. 


with adiabatic turn-on initial condition lim kF(7^ | 7^^(r^), r^) = 0. 


From ean. (IA23p and (|A2p the kinetic equation for the Ring DF F{IZ,t) is: 

1 


dF / 1 

^ + {^"iv 


[F, 4>(7^,r)] = 


N 


T(7^,7^'),4?(7^,7^^T) d7^^ 


(45) 


Fi2^(7^,7^^T) = kF(7^|7^^r)F(7^^T) + kF(7^'| 7^, r) F(7^, r) 

+ j W{IZ\IZ'',T)W{IZ'\IZ'',T)F{IZ'',T)dIZ''. (46) 

Eg ns. (145 D and (j46p are the basic equations governing the RR of a non-relativistic and isolated 
Keplerian stellar system. In the Appendix we discuss the 0(e) fluctuations and argue that they 
are bounded and hence eqns. (j45p and (I46p form a closed and consistent set of equations in 5-dim 
7^-space. An important feature of the evolution of F{IZ,t), as described by eans. (l45p and (|46p . is 
the absence of the mass ratio e in all the physical quantities, except in the time variable r = et. As 
pointed out by RT96 this implies the following property of mass-invariance: if the mass of every star 
in a Keplerian cluster is changed from m* to m*i, then the collisional evolution of the transformed 
cluster is identical to the original one, so long as the time variable r is replaced by the new time 
variable ti = (m*i/m*)T . 
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16 Sridhar & Touma 


5 RESONANT RELAXATION 

Here we provide a self-contained description of the equations governing RR, after generalizing 
to include orbit-averaged relativistic corrections and external perturbations; these will, of course, 
break the property of mass-invariance of time evolution mentioned above. 


5.1 Collisionless limit 


The Ring mean-field Hamiltonian H(TZ, r) can be obtained by supplementing the Ring mean- 
field potential $(7^, r) with general relativistic corrections and external perturbations. Relativistic 
corrections to the gravity of the MBH cause precession of both the apse and node of a Gaussian 
Ring. The relativistic (secular) Hamiltonian governing the 1 PN Schwarzschild precession of the 
apses and the 1.5 PN Lense-Thirring precession of the apses and nodes for a spinning MBH is 
given by: 


R-'(I,L,L,) 


+ 

'(/,l,l,). 


H^^\I,L) 

HL 

with 

„ 3{GM,)^M, 

^ M ' 


^^1'5(/,L,L,) 


with 

^ 2{GM,fxM, 

ci M • 

(47) 


where 0 ^ x ^ 1 is the spin parameter. Here we have assumed that MBH spin angular momentum 
points along the 2 ;-axis; factors (M,/M) = have been included in the constants Bi and Ri.s, 
because we measure time in terms of the slow time t = et. External perturbations due to nuclear 
density cusps and/or slowly moving distant masses can also be included. Let r) be the 

orbit-averaged potential and eA®^’'(r) be the acceleration of the MBH. These external sources will 
contribute a (scaled) orbit-averaged tidal potential. 


t) = $"^‘(7^,r) + X(7^) • Ar(r) 

dm 


where ^(^) = 9 -= centre-of-mass of a Gaussian Ring. 

' 2tt 


(48) 

(49) 


Then the general mean-field Ring Hamiltonian is: 

i^(7^,T) = <h(7^,r) + + 4>*'^(7^,T), (50) 

where we note that the mean-field potential is, as given in ean. HAlSp . determined self-consistently 
by the Ring DF: 

<1>(7^, r) = / F{'R!,T)'^{TZ,TZ')dTZ', Ring mean-field potential; 


where 


^(7^,7^') = -GM, 


dm dm' 1 


27r 27r \r — r'\ ’ 

The (mean-field) Ring equations of motion are: 


“bare” inter-Ring potential. 


(51) 

(52) 


I = \/GM,a = constant, 

dL _ dH dg _ dH 

dg ’ 


dL, 


dH 


dh dH 


dr dg dr dL dr dh dr dL^ 

Therefore, in the collisionless limit, the appropriate generalization of ean. (|37p for the F(TZ,t) is: 

dF 

+ IF,H\ = 0. 


(53) 
is: 

(54) 
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Stellar Dynamics around a Massive Black Hole II 17 
This is the general form of the Ring CBE, whose properties are discussed in more detail in Paper 1. 


5.2 Kinetic equation for Resonant Relaxation 


We now need to generalize ean. (l45l) for the Ring wake function, and (H6l) for the Ring DF, to account 
for the change in the mean-field Hamiltonian from d>(77, r) to H(IZ,t). The new wake equation 
can be determ ined frq r n the collisionless response of the stellar system, by extending the gedanken 
experiment of iGilbertI (| 19681 1 for point mass stars to Gaussian Rings: from the N Rings that are 
distributed in 5-dim Ring space according to the DF F (IZ, r), select one Ring and place it at the 
phase space location IZ'] this will induce a small perturbation [1/N)W at every location IZ in the 
5-dim 77-space. Then the net perturbed DF at any instant t' < t is: 


Fi(77,r') = - 


F{IZ,t') 5{IZ-IZ'{t')) W{IZ\IZ'{t'),t') 


+ 


(55) 


N N N 

where IZ'{t') is the location at time t', of the Ring which arrives at IZ' at time r. Since N ^ 1 the 
perturbation is effectively infinitesimal, so Fi must satisfy the linearized CBE discussed in Paper I. 
This then gives the general Ring wake equation we seek: 


dW 


+ [W{IZ\IZ'{t'),t'),H{IZ,t')] + [F(77,r'), d-'"(77,77'(r'),T')] 


= [<hP(77,77'(r'),r'),E(77,r')] , 


for t' < r, 


with adiabatic turn-on initial condition lim W{IZ\IZ '= 0. (56) 

r'^—oo 

Here is the gravitational potential due to the wake which was derived in ()A20h : 

$"(77,77',t') = j W{IZ"\IZ',T')^iIZ,IZ")dIZ" , Ring wake potential. (57) 

‘hP is the difference between the “bare” inter-Ring interaction potential and the mean-field potential, 
which was derived in (IA21I1 : 

<hP(77,77',r') = T(77,77') — <h(77, r'). Ring perturbing potential. (58) 

As earlier the Ring wake function satisfies the two identities 

J W{IZ \ IZ', t) dIZ = 0 , Zero mass in the wake of Ring IZ'] 

J W{IZ \ IZ', t)F{IZ', t) dIZ' = 0 , Zero net wake at IZ due to all the Rings. (59) 

When we compare eqn. (156h with the earlier eqn. (l45h for the Ring wake function, we see that the 
only change is in the second term on the left hand side where $ has been replaced by H. This is 
expected because the mean-field Ring orbits are now governed by H. However on the right hand 
side, the perturbing potential 4>P = T — remains unchanged, because this is the driver of the 
irreversible statistical evolution of the system, and cannot be affected by deterministic agencies like 
relativistic corrections and external gravitational fields. 

The new kinetic equation for the DF F{IZ,t) is obtained by replacing, on the left hand side of 
ean. (l46p . the quantity (4> — ^/N) by {H — ^/N), because the ^/N term corrects for the fact that 
only (N — 1) of the N stars exert forces on any given star. The right hand side, which represents 
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irreversible statistical evolution, remains unaltered. Therefore the general Ring kinetic equation is: 


OF 


+ 


F,H- 


Hn,T) 

N 


1 

N 


T(7^,7^'),4?(7^,7^^r) d7^', 


F^^J{n,n',T) = w{n\n',T)F(n',T) + w{n'\n,T)F{n,T) 


+ 




(60) 


( 2 ) 

By eliminating and manipulating the resulting expressions, we can also write the Ring kinetic 
equation in a form where the dissipative and fluctuating contributions to the Ring collision term 
are displayed explicitly: 


OF 


+ 


F,H- 


<h(7^,r) 


N 


= + C^'^^[F] , 


'flue r 


(61a) 


j F{n,T)w{n'\n,T)] dn', (6ib) 

c»-[ri = 1/ [nn,n') + w(n\n',r)].n'. (eic) 


The kinetic equation — either ean. (l60p or eans. (l61ap - (l61cp — applies to the collisional evolution 
of F whenever it satisfies 


F{n,T),H{n,T) 


^(^,t) 

N 


~ 0, 


being of order 1/N or smaller. 


(62) 


Therefore F{TZ,t) may be thought of as passing through a sequence of collisionless equilibria in a 
quasi-static manner. 


Boundary conditions: Eqn. (|56P and either eqn. (l60h or eqns. (l61ap - (l61cp need to be supplied with 
boundary conditions in F-space. There are two main cases of interest: 


1. “Lossless” stellar systems in which the MBH is not considered to be a sink of stars. Then the 
DF is normalized as f F(F, r) dTZ = 1. Subject to this normalization F{TZ, r) can take any positive 
value at any location in F-space; in other words, the domain of F is all of F-space. Hence the 
domain of the wake function W{TZ \ TZ', r) is all of TZ and TZ' spaces. 

2. Stellar systems in which stars that come too close to the MBH are either tidally shredded or 
swallowed whole. Since the system loses stars, we can only require that the normalization of the 
DF at some initial time tq, be 

J F{TZ,to) dTZ = 1, Initial normalization with loss of stars to MBH. (63) 

At later times r > tq, we will have J F{TZ,t) dTZ ^ 1. The simplest model of the loss is the 

assumption that a star is lost to the MBH when its pericentre distance is smaller than some fixed 
value ric, the loss-cone radius. When the loss-cone is empty, stars belonging to the cluster must 
necessarily have pericentre radii a(l — e) larger than r\c. Hence the domain of F{TZ,t) is restricted 
to regions of F-space in which I and L are large enough: 


lie < T where 

Lic(I) < L ^ I, where 


he = VGM.rie ; 


Lle{I) 


he 


2 - 



21 1/2 


(64) 


As I increases from its minimum value of Ac, 


the function L\e{I) monotonically increases from 
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its minimum value of Ijc and approaches its asymptotic value of • The Ring DF and wake 

functions satisfy the “empty loss-cone” (or “absorbing”) boundary conditions: that for all r, 

F[TZ,t) = 0 for I ^ lie and L ^ Lic(/), 

W{n I 7^^ r) = 0 for I, ^ lie and L ^ Lic(/), L' < Lic(/')- (65) 

Eon. (1561) and either ean. dbOp or (I61ap - (l61cp . together with suitable boundary conditions on the 
DF and wake, are the fundamental equations governing Resonant Relaxation. 


6 DISCUSSION 

It should come as no surprise that there is a structural similarity between the Ring kinetic equation, 
and the Poisson Bracket form of Gilbert’s equationj^This invariance of form is a natural consequence 
of averaging over a (fast) angle variable, a procedure that preserves canonical structure in phase 
space. Gilbert’s equations apply to a general stellar system in 6-dim phase space, whereas the RR 
equations apply to a Keplerian stellar system orbiting a MBH in a reduced 5-dim Ring space. What 
we have demonstrated is that the Ring equations can be got from the 6-dim Gilbert equations by 
replacing (a) “point mass star” with “Gaussian Ring”, (b) all 6-dim PBs with 4-dim PBs, and (c) the 
time variable t by the slow time r = et. Thus we have corresponding (scaled) potentials: ^{IZ,IZ') 
which is the “bare” inter-Ring potential between Rings IZ and IZ'] the Ring mean-field potential 
<h(7?., r) at location IZ; the potential {JZ,1Z' ,t) felt by a Ring at IZ due to the wake of Ring IZ'; 
and ^^{IZ,IZ\t) = ^{IZ,IZ') — <I>(7^,r) which is the perturbing potential felt by a Ring at IZ due 
to a Ring at IZ' . All these potentials are related to the DF or the wake through formulae that are 
analogous to those in Gilbert’s theory. 

Ring Wake as the driver of RR: The kinetic equation in the BBGKY form ean. (l60ll tells us 
that irreversible collisional evolution is driven by the irreducible part of the 2-particle correlation: 

4y(7^,7^^r) = W{IZ\IZ',t)F{IZ',t) + W(7^'| 7^, r) F(7^, r) 


+ / W{IZ\IZ",T)W{IZ'\IZ",T)F{IZ",T)dIZ" 


( 66 ) 


(o') 

This decomposition of FM means that the irreducible 2-Ring correlation at {IZ,IZ') gets three 
kinds of contributions from the wake function: (a) The wake of IZ' at IZ; (b) The wake of IZ at IZ'; 
(c) The product of the wake values at the points IZ and IZ' of a third Ring at IZ", summed over 
all locations IZ". All three contributions come with suitable T-weighting. The third term accounts 
for the contribution of collective effects (“gravitational polarization”) to the microscopic processes 
driving RR. This is the nature of the full theory at 0{1/N). 

Correlations build through Ring-Ring collisions: very early times r —?> — oo, the wake W ^ 0 


which implies that F- 


( 2 ) 


( 2 ') 

0. As the wake of every Ring builds over time, so does FI ' so we can 


think of the wake as the fundamental driver of collisional evolution. The Ring wake ean. (l56p is: 
dW 


dr' 


+ [W{IZ\IZ'{t'),t'),H{IZ,t')] + [F{IZ,t'),^'"{IZ,IZ'{t'),t')] 


= [4>P(7^,7^'(r'),r'),T(7^,r')] 


for t' ^ r. 


with adiabatic turn-on initial condition lim W{IZ\IZ'{t'),t') = 0, 


^ Compare the Ring wake eqn. iHg with eqn. C3; the Ring kinetic equation in the BBGKY form eqn. IMli with eqn. m ; or 
the Ring kinetic equation in fluctuation-dissipation form eans. l|61a|l — Il61c|l with eans. l|15a|l — Ill5c|l . 
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where 

7^^ t') = j W{n" I 7^^ t') 7^") d7^" , Ring wake potential; 

,t') = Ring perturbing potential. 

Note that is linear in W, whereas ‘hP is independent oiW. If the right hand side happened to 

be zero, i.e. , F{TZ,t')] = 0, then the wake equation would be a linear integral equation which 
is homogeneous in W. The solution that is compatible with the adiabatic turn-on initial condition 
is Vk = 0 for all time. Since there is no wake, the collision term in ean. (l61al) - (l61cp vanishes. Then 
the Ring DF satisfies the Ring CBE, and there is no relaxation. Therefore 

[ <i>P(7^, 7^'(r'), r'), F(7^, r') ] = [ 'k(7^, 7^'(r0) - <I>(7^, r'), F(7^, r') ] (67) 

is the “source term” for the Ring wake function, and for the entire RR process. In general this 
drives changes in both the magnitude and direction of angular momenta through apsidal and nodal 
resonances, and there is no strict separation into scalar-RR and vector-RR. 

Time evolution: The Ring kinetic equation is for RR what the Boltzmann equation is to the 
kinetic theory of gases. Just as the Boltzmann equation follows the DF over times that are longer 
than the duration of collisions between gas molecules, the Ring kinetic equation tracks changes in the 
DF of the system over times that are longer than Tgec which is the typical coherence time of secular 
stellar encounters. Over times several T^ec ^ NTgec, the collision term 

makes only a small contribution to the change in F. Then the system is effectively collisionless and 
is well-approximated by the Ring CBE eon. (1541) . In this phase the Ring DF can display signihcant 
variations over times Tsec- these could be due to (a) secular instabilities, or (b) secular collisionless 
relaxation (i.e. violent relaxation of Gaussian Rings), or (c) excitation of secular modes by an 
external perturber, and possible loss of stars to the MBH. In the absence of continuing forcing 
on the secular time scale by external sources, it is expected that the stellar system would settle 
in a stationary state dF/dr = 0 which implies that \F,H] = 0. By the secular Jeans theorem 
of Paper I, F must be a function of the isolating integrals of motion of H. These Ring DFs are 
necessarily linearly stable to perturbations. 

Over a longer time scale 3> Tsec the effect of the collision term can no 

longer be ignored. The Ring DF will evolve slowly due to the slow accrual of the relaxing effect of 
angular momentum exchanges between the Rings. We can think of F as passing slowly through a 
sequence of secular Jeans equilibria for which [F, FI — <I>/A^] is of smaller order than 1/N. Then 
eqns. (l6Iap - (j61cp imply that dF/dr ~ ~ 0(I/A^). Therefore the time over 

which F deforms (quasi-steadily) by order unity is the RR time scale Tres = NT^ec = {N/£)T’^ep ■ 
During this phase of evolution we can allow for the external tidal potential to also vary slowly 
— as might be the case when a galactic nuclear cusp adjusts adiabatically to the slowly varying 
Ring DF. All through the relaxation process. Rings exchange angular momentum, so a Ring that 
has lost enough angular momentum to have its periapse decrease to rjc will be lost to the MBH. 
It is of great interest to determine the mass, energy and angular momentum lost by the Keplerian 
star cluster to the MBH. 


7 CONCLUSIONS 

We have derived the fundamental kinetic equation governing the Resonant Relaxation (RR) of low 
mass stellar systems around a massive black hole (MBH). This includes the effects of stellar self¬ 
gravity, general relativistic corrections up to 1.5 post-Newtonian (PN) order, as well as external 
deterministic sources of gravity that vary on the secular time scale. We considered a stellar system 
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consisting of stars of equal mass, but this is not a serious limitation. Mass-segregation effects in RR 
can be studied by straightforward generalization to a system with a range of stellar masses. We have 
seen in the discussion of the previous section that the wake function is the driver of RR, and that 
the source term for the wake has contributions from both apsidal and nodal resonances which may 
be equally important in general geometries. Therefore the traditional, physically-motivated split of 
RR into scalar-RR and vector-RR dissolves, in general, into a seamless blend of both processes, 
with the later dominating when fast apse precession promotes apse-averaged dy namics. We noted 
the structural similarity between the Ring kinetic equation and Gilbert’s equation. Hevvaerts ( 201Cll ) 
proved an H-theorem for the latter, for an integrable stellar system admitting global action-angle 
variables. It seems reasonable to expect that a similar exercise would be successful for RR. Our 
theory of RR is valid for stellar systems with arbitrary geometry and figure rotation, so long as 
these are dynamically stable quasi-steady secular equilibria. The orbital structure can be regular, 
chaotic, or mixedjl the formalism presented in this paper applies to all of these. 

Below is a brief guide for the application of the Ring kinetic equation to any Keplerian stellar 
system. These are followed by remarks on future directions and useful ways of comparison with 
numerical simulations. 

1. Secure an “explicit” form for the orbit-averaged Poisson kernel 'k in eqn. (f52ll . Although desir¬ 
able, a closed form expression is difficult to come by; more often than not the orbit averaged “bare” 
Ring potential is expressed as a Fourier series in the apses and nodes, with the Fourier coefficients 
depending on the eccentricities and the inclinations (or related Delaunay variables). One could also 
work with a pre-computed table of the kernel in the physical limit of interest. 

2. Using 'k construct an equilibrium pair, {F , , of DF F and mean-field potential, <k, which 

is related to F through the integral of eqn. (151 jl . Secular dynamics is governed by H, the mean-field 
Hamiltonian of eqn. flhOj) . which is the sum of the mean-field potential <k, and contributions from 
1.5 PN general relativistic corrections and any external gravitational potential. By the secular Jeans 
theorem of Paper f, F must be a function of H. In general one needs to solve an integral equation; it 
would be very interesting to see the development of Schwarzschild-like iterative numerical methods. 
Moreover it is also necessary that the F{H) be dynamically stable. 

3. With orbit-averaged kernel in hand, and an initially stable equilibrium model, solve the wake 
equation (I56p for W, with adiabatic turn on initial condition. 

4. Use the pair {F , VF} to compute the collision integral, which is the right hand side of the Ring 
kinetic equation (l61aD - (l61ch driving relaxation. 

5. Solve the kinetic equation (Eq. I61ap - (l61cp for F ■. the global solution of eqn. (l61ap couples all 
the three steps given above, with the Ring mean-field potential $ recovered from E, then wake W 
updated with F and <k as one progresses towards a relaxed stated 

The steps above are explicitly laid out in a companion paper [^ridhar &: Touma ( 2016: Paper fil l, 
where we formulate and s tudy the physical kinetics of the R R of zero-thickness, flat, axisymmetric 
discs. Working with the IPolvachenko fc ShukhmanI (Il982l ) approximation of neglecting “gravita¬ 
tional polarization” we derive a Fokker-Planck equation for the DF with diffusion coefficients that 
are given self-consistently in terms of contributions from apsidal resonances between pairs of stellar 


® Averaging over the fast Kepler orbital phase is independent of whether the resulting secular dynamics is integrable or not; 
this is as true for planetary systems as Keplerian stellar systems. When the secular dynamics is non-integrable, chaotic diffusion 
through broken tori is restricted to the 4—dim {L, Lz, 5 -, h} subspace of the full 6 —dim phase space. 

^ In the course of rela.xation, the linear secular stability of the evolving quasi—equilibria must be monitored. In case one 
encounters a dynamical instability, fast (i.e. on the secular time scale) collisionless processes will intervene. Once the instability 
has saturated and the system settled in a new secular collisionless equilibrium, the subsequent RR evolution can be resumed 
by following steps 2-5. 
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orbits. For ‘lossless’ stellar discs we prove an H-theorem and conservation laws for the disc mass, 
energy and angular momentum. When stars are lost to the MBH through an empty loss-cone, the 
loss rates of mass, energy and angular momentum can be expressed in terms of the diffusion flux 
at the loss-cone boundary. 

It is hard to overemphasize the need for comparison with numerical simulations, and we can 
see three different routes of exploration: (a) The diffusion coefficients and angular momentum 
diffusion statistics predicted by our RR theory could be compar ed with results obta ined from 
refined N-body simulations designed for that same purpose, such as iMerritt et al.l ((20111) for RR in 
a spherical cluster; (b) The Ring kinetic equation (|61ap - (|61c|) can be solved for the thermally relaxed 
equilibrium of any chosen initial stable equilibrium. The equilibrium and associated statistics could 
be compared with the N-body evolution of (a sample of) the same initial state; (c) The time 
evolving DF of our formalism can be compared with related solutions of various Fokker-Pl anck 
equations proposed by others ( Hopman fc Alexander 2006 : Madigan, Hopman fc LevinI 2011 1. ad 
hoc though their assumptions may be. This should allow calibration of the performance of ad hoc 
recipes for RR vis-a-vis the fundamental formulation we have presented. 
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APPENDIX A: DERIVATION OF THE ORBIT-AVERAGED KINETIC 
EQUATION 


Here we derive the orbit-averaged kinetic equation, beginning with the exact Gilbert equations 
eqns. (l3T]l and (I32p given in Delaunay variables. Orbit-averaging is a systematic procedure, of per¬ 
turbative development of these equations in the small parameter e. 

The first step is to substitute in (I32p the multiple-scale forms, eqn. (j43)) for the DF / and eqn. (j44ll 
for the wake W. Then all the 0(1) terms vanish. The 0(e) terms give the following equation for 
F{IZ,t) and all the /„(77.,r): 


OF 


+ ^ inOk fn exp [inw] + 
n^O 




ip + r-A, 


(6) 


where 

4?(7^,7^^r) 


1 

N 


p{v,v'),F^^>{n,n',T) 


(6) 


d7^ 


,dw' 


VF(7^|7^',r)F(7^',r) + 1F(7^'| 7^, r) F(7^, r) 


(Al) 


+ y lT(7^|7^",r)lT(7^'|7^",r)F(7^",T)d7^", (A2) 

is the irreducible part of the 2-Ring correlation function. Note that only the 0(1) part of the wake 
W contributes directly to the 0(e) kinetic theory. The Poisson kernel = —GM,/|r — r'\ 

is a specified 0(1) function of its arguments and is independent of e. The other quantities, (p and 
A,, are needed only to 0(1) accuracy. From ean. (|28cp we have to 0(1): 

A.(r) = GM. J F{n, r) d7^ ^ ^ ^ = 0 , (A3) 

because f dwr/r^ = 0, by the conservation of angular momentum along a Kepler orbit. Recalling 
that X = eA, we see that, as in the collisionless theory of Paper I, the MBH does not accelerate to 
0(e) accuracy. From eqn. (l28bp the rescaled mean potential to 0(1) is: 

(^(P,r) = j F{n\T)dn' (A4) 
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The PB on the right hand side of ean. fjAip are over the T> variables, so the integral over w' can 
be taken inside the PB. Since p{T),T)') is the only quantity inside the PB that depends on w', the 
tc^-integral operates only on p{'D,'D'). Then ean. HAip simplihes to: 

— + ^ in^k/n exp [inu;] + f 1 - — j [F, r) ](g) = 


1 

iv 


/ 



dW . 


Averaging over w, we get the following equation for F{TZ,t): 




[F,<p{V,t)] 


dw 

(6)^ = 


(A5) 


1 

N 



p{V,V') 


dm dm' 
27r 27r 


Fi?(P,P',r) 



(A6) 


Subtracting ean. llAbp from (IA5I) . and dropping the extremely small terms proportional to 1/N in 
comparison to 0(1) terms, we get the following equation for the /„(P, r): 

[F, y?(P,r)](g) ^ = 0, (A7) 

nj=v 

with solution 

/n(F,r) = [F, (^(P,r)](g)exp[-inm] n ^ 0 . (A8) 

Hence all the /„(F, r) for n^O are slaved to the quasi-steady collisional evolution of F(F, r). As 
in Paper I, the function /o(F, r) remains undetermined in our 0(e) calculations. 

Similarly, when the expansions of eans. (l43p and (14411 for the DF and the wake are substituted in 
eqn. m, all the 0(1) terms vanish. The 0(e) terms give the following equation for W{TZ \ TZ'{t'),t') 
and all the ci7n(F. | P'(t'), r'): 


dW 

W 


+ ^ inHk exp [inm] + [ IT , (^(P, r') ] ^g^ 
n^O 


+ [F(P,r'),^-(P,P'(r'),r')](g) 


= , F{n,T')]^^^ , forr' ^ r. 


(A9) 


Note that we have written the rescaled wake potential as (/^"(P,P',r), instead of (^"(P,P',t). 
This is because, when the expansion of eon. (1441) is substituted in ean. (l28dl) . we find that (/?* is 
independent of m' to 0(1) accuracy: 

p^{v,n',T) = j w{n"\n',T)dn" ^ p{v,v")^ . (aio) 

In ean. (IA9p P'(t') = {P'(r'), m'(r')} is the location of the star in 6-dim phase space at time r', 
which arrives at the phase space location P' = {TZ', m'} at time r. In general this orbit is governed 
by the Hamiltonian Forg of ean. (l2^ . but we need it in ean. (IA9p only to 0(1) accuracy. This is 
just the secular collisionless dynamics described above by eqns. fjdTI) and (I42p : TZ'{t') is such that 
I' = constant and {F'(r'), L^(r'); g'{ t'), /i'(r')} obey the Hamiltonian equations of motion eqn. (l4ip 
with 4>(P',r') acting as the Hamiltonian. Integrating eon. (1421) over time, we have 


m'(r') = m' + 


Hk(F) 


(r'-r) + ^ 


54>(P'(t"),t" 

m 


dr" 


(All) 


The Ring orbit TZ'{t') is independent of which is slaved to it. 

Since only the IT part of the wake contributes to the equation for F and /„, we can obtain 
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an equation for it by averaging ean. (IA9|l over both w and w'. All the terms with Wn vanish in the 
re-averaging, and play no further role. Then the equation for WiJZ\IZ' is: 


dW 


+ 




dre dw' 
27r 27r 


r /X /X /X n dre dre' 

[¥pP(P,P(T),r), F(7^,T)](g) 


for t' < 


T. 


with adiabatic turn-on initial condition lim W{IZ\IZ '= 0. 

r'^—oo 


(A12) 


Equation (IA12P for W, ean. (IA6h for F, ean. ()A2p for and ean. ()A8h for the /„ are the orbit- 
averaged equations describing the 0(e) kinetic theory. However they are in a raw form and need to 
be processed further such that all direct dependences on the orbital phases w and w' are eliminated. 
In other words, we seek a kinetic description that is self-consistent and closed in 5-dim 7^-space. 
This is done in below by (a) expressing the 6-dim PBs in terms of 4-dim PBs, and (b) defining 
suitable potential functions in 7^-space. 

For any two phase space functions, Xi(^) X 2 {D), the 6-dim PB can be expressed in terms 

of the 4-dim PB (without subscript) as 


[Xl , X2](6) - 


_ dxi dx 2 
dw dl dl dw 


+ [Xi , X 2 ] 


(A13) 


For functions Xi{H) and X 2 iH) , the 6-dim PB equals the 4-dim PB: [xi, X2](6) = [Xi > X 2 ]- 
We begin with the 3 PBs in the wake eqn. (IA12h : 


[y,P(V,V\T'),T'),Fin,T')] 



dW dip 

(6) 

dl dw 


dF dip"** 

(6) “ 

dl dw 


, dp>^ dW 

(6) 

' dw dl 


+ [IT(7^|7^'(r'),r'), 

+ [F(7^,r'), ¥^"'(P,7^'(T'),r')] 

+ [ipP(V,V'(T'),T'),F(n,T')] 


The first term on the right hand side of all three equations can be written in the form 9{ }/dw, and 
will vanish when averaged over w. So the three double-integrals over the 6-dim PBs in ean. (IA12h 
reduce to double-integrals over the 4-dim PBs. Since the 4-dim PBs act only in 7^-space, the 
integrals over w and w' can be taken inside the 4-dim PBs. Of the two functions in each of the PBs, 
it is only the potentials ip, ip"*' and ip^ that depend on the orbital phases w and w'. In particular ip 
and ip"** are functions of w but independent of w', so we get: 


[w(n\n'(T'),T'),ip(v,T')] 


drc dtc' 

(6) 


w{n\n'(T')y), 


[F(K,T').v"(P.K'(T').r')](e, 



However, ip^ depends on both w and w', so we are left with a double integral inside the 4-dim PB: 


[ip^(V,V'{T'),T'),F{n,T')] 


dm dm' 
(6) 27r 27r 


jj^P(V,V\T'),T')^^,F(n,T') 


In the integral over m' on the right hand side, we note that m' occurs in the integrand only through 
the function m'(r') which is given in ean. dAlip . Since m' appears here only as a linear additive 
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term, we can replace w'{t') by w' in the potential 




dw dw' 


r,/„ / /.dtcdtc' ,, 

ip^{V,n (r ),'w ,T ) —— , F{n,T) 


2tt 2tt 


' (®) 27r 27r 

where the dependence of on the arguments TZ'{t') and w' is now indicated explicitly. 
Then eqn. (jA12p for W{TZ\TZ'{t'),t') can be written as: 

dW 


dr' 


+ 




ipP{V,n'{T'),w',T')^^ , T(7^,r') 


for t' ^ r. 


27r 27r 

with adiabatic turn-on initial condition lim W{TZ\TZ'= 0. 


(A14) 


We now turn to ean. (lA6p for the Ring DF F{TZ,t). Since 

dF dip 
dl dw 

we have 


f)cn 

[F(7^,r), (/::(P,r)](g) = + [T(7^, r), r) ] , 


Also 


j> [F{n,T), ip{V,T)]^^^^^ = j>[F{n,T),ip{V,T)]^ = F(7^,r), 


^i^dwdw' ^(2),^ ^ 

p{v,v')——,Fg\n,n\T) 


Then ean. ()A6p for F{TZ,t) can be written as: 


dF _ - 

^ + b“jv 


F, i>v{'D.r)^ 


( 6 ) 

1 

N 


j)j> p{V,V' 


dm dw' ^( 2 ) 
27r 27r ’ ■ i"" 


, Fl^>in,n',T) 


Similarly, the solution eqn. flASp for the 0(e) correction to the DF can be simplified to: 
1 dF f , dw 


d7^^ 

(A15) 


/n(7^,T) = 


Dk dl 


j) (p(T>, t) exp [—intc] 


27r 


+ 


nDk 


F , (f ip('D, t) exp [—inrc] — 

J 27r_ 


n 7 ^ 0. 


(A16) 


We have now eliminated all the 6-dim PBs and written ean. (IA14l) for VF, ean. HAlSh for F and 
ean. (IA16p for the /„ in terms of 4-dim PBs operating in 7^-space. Only the integrals of various 
potential functions over w and w' remain to be simplified. Using eans. (IA4ll for ip , (jAlOh for (/?* and 
ipP = p — ip ^ the integrals needed are: 

, dtc dw' 


2ti 2ti 


<i>(7^,r) = j^ip(V,T)^ = j F(n',T)dn' j^j^piV,V' 
iPniF, t) j) ifiV, t) exp [—intc] ^ = J F{TZ', r) dTZ' j)j) p(F, V') exp [—\nw] 

4•"(7^,7^^T) = j) <p'''{v,n',T)^ = J w(n”\n',T)dn" p(v,v") 


dw dw' 
271 27r ’ 

dm dm" 


27r 27r ’ 




dm dm" 
27r 27r 


dm 




A quantity that is common to all these integrals is the partially orbit-averaged Poisson kernel. 
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whose Fourier series we write as: 


_ .. dw' 

II 

+ y^'!/’n(77,77')exp[: 


n^O 

T(77,77') 

f 

f dtc dw' 1 

= -GM, (h 

/ 27r 27r r — r'\ ’ 

'0n(77,77') 

I 

II 

/ dtc dw' exp [—inw] 
J 27r 27r r — r'\ 


mw\ 


n / 0, 


(A17a) 

(A17b) 

(A17c) 


are the Ring-Ring interaction potential functions introduced in Paper I. The functions T and ipn 
will be treated as known functions of their arguments. Using eans. (IA17ap - (IA17cp . we get 


<h(77,r) 

= j F{n',T 

)T(77,77')d77', 

(A18) 


= j F{n',T 

)V^„(77,77')d77', 

(A19) 

4>"(77,77',r) 

= J W{IZ" 

77',r)T(77,77")d77", 

(A20) 

4>p(77,77',t) 

= T(77,77') 

- ‘^>(77,r). 

(A21) 


We are now ready to cast ean. ()A14ll for W, eans. ()A15|l for F and ean. ()A16p in a form that is 
determined entirely by 77-space quantities. The equation for the wake function W (77 | is: 


dW 


+ [W(77|77'(r'),r'),<h(77,T')] + [ F(77, r'), ^"(77,77'(r'), r') ] 


= [ 4>P(77,77'(r'), r'), 7^(77,/)] , for r' ^ r. 


with adiabatic turn-on initial condition lim kF(77 | 77'(r'), r') = 0. 


(A22) 


The kinetic equation for the Ring DF F{IZ,t) is: 


dF 


+ 1 - 


[F, ‘h(77,r)] = 


1 

N 


T(77,77'),7^i?(77,77',r) 


d77'. 


(A23) 


The Fourier coefficients of the 0(e) fluctuations in the DF fn are given by: 

1 f dF i 1 

/n(77, t) = n ^ ^ J ’ n ^ 0 . (A24) 

As earlier the function /o(77, r) remains undetermined by our 0(e) theory. Using eqn. (IA24jl for the 

fn in ean. (l43l) . we get the following expression for the full DF to first order in e: 

fi'F>,r) = ^F(77,t) + ^ [F, expPinu;] + ^/o(77,t). (A25) 


Once the Ring DF F(77, r) has been determined by solving eqns. (|45p and (I46I) . the 0(e) fluctuations 
in w are completely determined by eqn. (IA25p . To determine the function /o(77, r) it is necessary to 
work to higher order in e. With this caveat in mind, the orbit-averaged kinetic description presented 
here is self-consistent and closed in 5-dim 77-space. 
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